function j=dfun(x) 
F=fun();
jacobi=[diff( F,'X') diff( F,'Y') diff( F,'Z') diff( F,'A') diff( F,'B') diff( F,'G')]
j=jacobi

